





c^ 



o 



X 



Hole Pockets in the Doped 2D Hubbard Model 

Adriana Moreo and Daniel Duffy 
Department of Physics, National High Magnetic Field Lab and MARTECH, Florida State 

University, Tallahassee, FL 32306, USA 
(February 6, 2008) 

Abstract 



^ . Tlie electronic momentum distribution n(k) of the two dimensional Hub- 

^L) • bard model is studied for different values of the coupling U/t, electronic den- 

Q . sity (n), and temperature, using quantum Monte Carlo techniques. A detailed 

On ■ analysis of the data on 8 x 8 clusters shows that features consistent with hole 



pockets at momenta k = (ztS, it^) appear as the system is doped away from 



1-^ ■ half-filling. Our results are consistent with recent experimental data for the 

O ■ cuprates discussed by Aebi et al. (Phys. Rev. Lett. 72, 2757 (1994)). In the 



range of couplings studied, the depth of the pockets is maximum at (n) ~ 0.9, 
and it increases with decreasing temperature. The apparent absence of hole 
pockets in previous numerical studies of this model is explained. 



Typeset using REVT^ 



One of the simplest models often used to describe the physics of the high temperature 
cuprate superconductors is defined by the two dimensional (2D) Hubbard Hamiltonian. Sev- 
eral normal state properties of these materials, such as antiferromagnetism, the behavior of 
the magnetic susceptibility with temperature and density, the optical conductivity and oth- 
ers, are qualitatively reproduced by this model. [|I| Although superconductivity has not been 
observed in Monte Carlo simulations p, self-consistent techniques suggest its existence at 
low temperatures in the dx2_y2 channel. P] In order to understand the pairing mechanism 
leading to such superconducting instability it is important to study the shape of the Fermi 
surface in this model, and compare it against angle-resolved photoemission experiments for 
the cuprates. Theories based on antiferromagnetism (AF) explain in a natural way the 
apparent d-wave symmetry of the condensate, but they also suggest the existence of "hole 
pockets" in n(k). [|] Since for some time evidence of a large electron-like Fermi surface has 
been reported in experiments 0, ideas based on AF have been questioned. Moreover, in 
agreement with these experiments, theoretical studies of the Hubbard ||^ and t-J Hamilto- 
nians 0, reported electron-like Fermi surfaces very similar in shape to the non- interacting 
case, even at low hole densities where antiferromagnetic correlations should still be strong 
in the system. Possible many-body effects have been invoked to explain such a behavior. 

However, on the experimental side, these results have been recently challenged by Aebi et 
al. PI Using a photoemission technique that allows the mapping of the whole Fermi surface, 
these authors have shown the existence of hole pockets in their data for Bi2212 consistent 
with a rigid band filling of the hole states predicted by antiferromagnetically based mean- 
field approximations. [Q In parallel to the experimental results, recent theoretical arguments 
by Dagotto et al. ^ conjectured that the large Fermi surface of the 2D Hubbard model 
observed at finite temperature can be made compatible with the existence of hole pockets. 
Their reasoning is that unless the temperature of the simulation is smaller than the energy 
difference of holes states at momenta k = (±f , ±f ) and (0, vr), (vr, 0) (which they found to be 
very small at U/t = 10), the hole pocket effect is washed out for thermodynamical reasons. 
These authors suggested that studies at stronger couplings, lower temperatures and larger 



lattices should provide evidence for the existence of hole pockets. 

Motivated by these challenging experimental and theoretical ideas, in this paper a careful 
analysis of Quantum Monte Carlo data for the 2D Hubbard model is discussed. We conclude 
that numerical evidence indeed supports the existence of hole pockets at k = (±f , ±f ) at 
low hole density. Thus, our results reopen the possibility that the normal state properties 
of the high temperature superconductors at half-filling can be qualitatively approximated 
by a band filling of the states obtained in the mean-field spin- density- wave approximation. 
The Hubbard Hamiltonian is given by 

H = -t E (cU.- + h.c.)+UE(niT-l/2)(nu-l/2)+/.Eni,^> (1) 

<ij>,cr i i,(T 

where C;^ creates an electron at site i with spin projection a^ Uio- is the number operator, 
and the sum (ij) runs over pairs of nearest neighbor lattice sites. U is the on site coulombic 
repulsion, t the hopping amplitude, and /i is the chemical potential. This Hamiltonian will 
be analyzed using standard Quantum Monte Carlo techniques. |[10| Clusters with up to 8 x 8 
sites will be studied varying the temperature T, filling (n), and coupling U/t. The notorious 
sign-problem prevents us from working at very low temperatures and large couplings. Those 
used in this paper correspond to the lowest T and largest U/t that can be reliably reached 



with confidence using the QMC algorithm. ||Tl| We will show below that they are enough to 
observe hole pockets in the Hubbard model. 

The first issue to analyze is whether antiferromagnetic correlations are strong in the 
system away from half-filling, since this is a necessary condition for the existence of hole 
pockets induced by AF. In Fig.l-a, the spin-spin correlation function C(r) = (— l)l^l(SiSi_^j.) 
is presented for U/t = 6 on an 8 x 8 lattice at T = t/4. At half-filling, (n) = 1, (filled cir- 
cles) the correlations are strong and remain finite even at the largest possible distance clearly 
indicating the expected antiferromagnetic long range order (LRO) in the system. More in- 
teresting for our analysis is that at density (n) = 0.9 (open circles) where antiferromagnetic 
LRO does not exist, it is clear that the spin correlations are strong up to a distance of 
about 3 to 4 lattice spacings. Increasing further the hole density to (n) = 0.75 (triangles). 



antiferromagnetisin is further suppressed but it remains strong between nearest-neighbors 
spin showing that the moments are still well-formed even at this fairly large hole doping. 
When the coupling is reduced to the value that has been more widely analyzed by QMC 
techniques namely U/t = 4, a qualitatively similar behavior is found at the same temper- 
ature. The spin correlations at this coupling (shown in Fig.l-b) are less developed, but 
LRO is still observed at half-filling. Working at (n) = 0.9 short range antiferromagnetism is 
also robust as it occurs at higher couplings. Decreasing the temperature to T = t/6 (shown 
in Fig.l-c) the spin correlations at U/t = 4 become as strong as for U/t = 6 at T = t/4. 
Then, Figs.l-a,b,c have shown that at a realistic density of (n) = 0.9, the antiferromagnetic 
correlations are sufficiently developed that they may originate hole-pockets in the Fermi 
surface. Eventually, these pockets should evolve into a noninteracting-like Fermi surface as 
the electronic density is further decreased from half-filling. 

Since the antiferromagnetic correlations are robust, then the natural question to discuss 
is why hole pockets have not been observed in previous studies (which were only carried out 
at U/t = 4). To address this question let us concentrate on the analysis of the momentum 
distribution n(k) = Y^a c^o-Ck.o-, where c^^ = Yjj s^'^'-'cj ,^. Note that n(k) has been previously 
measured |^0 but in these and some other papers, the shape of the Fermi surface was 
determined by calculating the position of the momenta where n(k) = 0.5. In principle, this 
procedure is correct. However, using such a convention here we show that the pockets could 
be overlooked due to finite temperature and/or small lattice effects. To understand this 
point, in Fig.2-a n(k) versus k is shown along the kx = ky direction for U/t = 6 at T = t/4 
on an 8 X 8 cluster and density (n) = 0.9 (open squares). Comparing the interacting results 
with noninteracting U/t = data for the same parameters (filled squares) we observe that, 
in both cases, n(k) becomes 0.5 at approximately the same momentum. Repeating the 
calculation along other directions in momentum space, and if the n(k) k, 0.5 criterion 
is still used, it is thus understandable why a Fermi surface resembling a non-interacting 
system was obtained in previous studies. However, this does not rule out the existence of 
pockets. The crosses shown in Fig.2-a correspond to spin density wave mean-field (SDW- 



MF) results. [Q In this approximation, an antiferromagnetic state is used which effectively 
produces a 2 x 2 unit cell. The mean-field Hamiltonian is diagonalized producing conduction 
and valence bands separated by the antiferromagnetic gap. The energy levels are given by 



Ek = iy Cfc + A^, where ek = —2t{coskx + cosky), and A is found using a self-consistent 
equation. At half-filling, the valence band is filled, and n(k) = |(1 — ^) is in very good 



agreement with the numerical data |T2| as it is shown in Fig.2-b. Now, let us assume that a 
robust antiferromagnetism survives the introduction of hole doping, as it was suggested by 
the results of Fig.l-a,b,c. Quasiparticles are removed from the top of the valence band to 
mimic the presence of doping. n(k) now is given by, 

n(k) = -(l-— )(— ^-^^— ), (2) 

where /i is selected such that the density is (n). In Fig.2a-b the SDW mean-field results are 
compared against the Monte Carlo data. The agreement along the kx = ky direction both 
at half-filling and at finite hole density is remarkable. We have explicitly checked that only 
at densities where the antiferromagnetic correlations become of one lattice spacing or less 
((n) ^ 0.75), the comparison between the numerical data and the SDW results deteriorates. 
Reducing further the density, the QMC data converges smoothly to a weakly interacting gas 
of electrons. 

The agreement between a theory based on strong AF correlations and the numerical 
data is consistent with the strong AF correlations shown in Fig.l. However, the main 
issue addressed in this paper remains paradoxical, i.e. if the numerical data and the SDW 
approximation are in good agreement, why there is no hole pocket in the occupation number 
mean value of Fig.2-a? To understand this point, it is instructive to study the effects that 
high temperatures and a finite lattice produce in n(k) in the mean-field approximation. In 
Fig.3-a we show the mean field n(k) along the k = kx = ky direction on a 20 x 20 lattice, for 
U/t = 6 and (n) = 0.9. At T = there is a clear pocket-like feature at {^, |), very different 
from that observed numerically in Fig. 2a. However, as the temperature increases the pocket 



becomes less pronounced and for temperatures T = t/6 and t/4, which are the values most 
commonly studied with Monte Carlo techniques, it has all but disappeared. In Fig.2-a the 
crosses are the mean-field results and it is clear that they are in much better agreement with 
the Monte Carlo data than the U/t = points. This shows that the existence of pockets 
along the k = kx = ky direction is a feature very difficult to see and quickly washed out 
by high temperature effects [Q due to the rapid change in n(k) where the pockets are 
expected. In Fig.3-b, we show the region on the U/t-T/t plane where a pocket, defined as a 
local minimum in the momentum distribution, is expected at (|, ^) according to the SDW 
mean-field results at (n) = 0.9. In the same plot the crosses indicate the lowest temperatures 
that can be presently reached for each U/t value with the QMC algorithm, the limitation 
being caused mainly by the well-known sign-problem. Clearly, all the crosses lie in a region 
where pockets are not observed in the mean-field approximation. Thus, the absence of a 
local minimum in our results of Fig.2-a are not incompatible with hole pockets, and it is 
caused by finite temperature and small lattice effects. 

The most striking evidence showing that the interacting system is different from the 
non-interacting one is obtained analyzing the line from k = (0, vr) to (tt, 0) (i.e. from Y to 
X). The SDW-MF predicts a constant momentum distribution along cosk^ + cosky = even 
at finite coupling, but this is a spurious degeneracy of this simple approximation and thus 
it is not useful to contrast against the QMC data. The noninteracting system U/t = 0, also 
has n(k) constant along this line at any density and temperature, and it is given by 

"^^^^ = e/3(^k-M) + 1 = e-f^'' + 1 ■ ^^^ 

At half-filling /i = and n(k) along this line is 0.5. This symmetry allows us to study pockets 
in the interacting system by looking for a reduction along the X-Y line of the momentum 
distribution as the coupling is increased. If n(|, |) is the smallest along this line, it would 
be an evidence in favor of hole pockets. In Fig.4-a n(k) vs k along the (tt, 0) — (0, tt) line 
is presented, using an 8 x 8 cluster, U/t = 6, and T = t/4. In the half-filled case (filled 
circles) n(k) is constant even for an interacting system. However, at density (n) = 0.9 
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(open circles) there is a clear minimum (pocket) at k = (|, |). The minimum lasts as long 
as antiferromagnetism is present in the system, i.e. at density (n) < 0.75 or smaller the 
minima are no longer observed giving support to our interpretation of this feature as hole 
pockets caused by AF correlations. 

In Fig.4-b, similar results are presented for U/t = 4. Here the pockets are smaller than 
at U/t = 6 because the AF order is not so strong but the qualitative behavior is the same. 
The maximum depth occurs again for (n) = 0.9. In Fig.4-c we show the pocket at several 
couplings U/t ranging from to 6 at a fixed density (n) = 0.9. As the coupling increases the 
pocket becomes deeper which is intuitively understandable using the AF picture. Finally, 
in Fig.4-d the dependence with temperature at U/t = 4 and (n) = 0.9 is presented. As 
with the coupling dependence, the pocket becomes deeper as the temperature is lowered 
from T = t/2 to t/6. ||14[ It is clear that the behavior of the pocket is correlated with the 
behavior of the magnetic correlations presented in Fig.l. In other words, when the spin 
correlations are strong, the pockets appear in the spectra. When the correlations increase, 
due to an increase in coupling or decrease in temperature, the pockets become deeper. But 
when the AF correlations are of the order of one lattice spacing, as for (n) < 0.75, the 
pockets are no longer observed, and n(k) becomes constant along the (vr, 0), (0, vr) direction 
as in the noninteracting case. 

In this paper, we have presented evidence that in the 2D Hubbard model the strong 
antiferromagnetic fluctuations produce hole pockets at momenta k = (±|, ±|). Monitoring 
the behavior of n(k) along the line (vr, 0), (0, it), where the electronic momentum distribution 
is constant in the noninteracting system, we have observed a minimum at (f , f ) compatible 
with the existence of a hole pocket. Along the kx = ky direction no local minimum is 
numerically observed but by comparing our results to a SDW-MF approximation we have 
shown that this feature is due to finite temperature effects, as predicted by Dagotto et al. 
For the values of U/t here studied, i.e. 4 and 6, the pockets exist for {n) > 0.75 and the 
maximum effect is observed at (n) = 0.9. Larger couplings and lower temperatures would 
certainly increase the range of existence of these hole pockets. Thus, we have shown that the 



recent experiments of Aebi et al. are compatible with results obtained for the 2D Hubbard 
model giving support to theoretical descriptions of the cuprates based on antiferromagnetic 
correlations. 
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Figure Captions 

1. (a) Spin-spin correlation C(r) = (SJSJ+r)!-!)'""' for U/t = 6, T = t/4 on an 8x8 lattice 
at different fillings. The error bars are of the size of the dots. 

(b) Same as (a) but for U/t = 4. 

(c) Same as (b) but for T = t/6. 

2. (a) n(k) as a function of momentum along the kx = ky diagonal direction on an 8 x 8 
lattice for U/t = 6, T = t/4 at (n) = 0.9. The open squares are Monte Carlo results, 
the filled squares are U/t = results, and the crosses correspond to the SDW mean 
field approximation (see text). 

(b) Same as (a) for (n) = 1.0 (half-filling). 

3. (a) Mean field values of n(k) vs momenta along the diagonal kx = ky direction using 
a 20 X 20 lattice, U/t = 6, (n) = 0.9 for T/t = 0, 1/50, 1/20, 1/6 and 1/4. 

(b) Mean field determination of the region where pockets, define as a relative minimum 
along the kx = ky direction, are observed in the U/t vs T/t plane for (n) = 0.9 on a 
20 X 20 cluster. The crosses indicate the lowest temperatures that can be reached 
with our Monte Carlo technique for different values of the coupling, due to the sign 
problem. 
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4. (a) Momentum distribution n(k) as a function of momenta along the (tt, 0) to (0, vr) 
direction on an 8 x 8 cluster for U/t = 6, T = t/4 at (n) = 1 (filled circles), (n) = 0.9 
(open circles) and (n) =0.75 (filled triangles). 

(b) Same as (a) for U/t = 4. 

(c) Momentum distribution n(k) as a function of momentum along the (vr, 0) to (0, tt) 
direction on an 8 x 8 lattice at (n) = 0.9 and T/t = 1/4 for (from above to below) 
U/t = 6,4,2 andO. 

(d) Momentum distribution n(k) as a function of momentum along the (tt, 0) to (0, tt) 
direction on an 8 x 8 lattice at (n) = 0.9 and U/t = 4 for temperatures (from above 
to below) T = t/2, t/3, t/4 and t/6. 
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